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Gamma-ray variability of radio-loud narrow-line Seyfert 1 
galaxies 
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ABSTRACT 

The recent detection of 7-ray emission from four radio-loud narrow-line Seyfert 1 
galaxies suggests that the engine driving the AGN activity of these objects share some 
similarities with that of blazars, namely the presence of a 7-ray emitting, variable, jet of 
plasma closely aligned to the line of sight. In this work we analyze the 7-ray light curves 
of the four radio-loud narrow-line Seyfert 1 galaxies for which high-energy gamma-ray 
emission has been discovered by Fermi/h AT, in order to study their variability. We 
find significant flux variability in all the sources. This allows us to exclude a starburst 
origin of the 7-ray photons and confirms the presence of a relativistic jet. Furthermore 
we estimate the minimum e-folding variability timescale (3 - 30 days) and infer an 
upper limit for the size of the emitting region (0.2 - 2 pc, assuming a relativistic 
Doppler factor S — 10 and a jet aperture of 6 — 0.1 rad). 

Key words: galaxies: jets - galaxies: Seyfert - galaxies: individual: PMN J0948+0022 
- 1H 0323+342 - PKS 1502+036 - PKS 2004-447 - gamma-rays: observations. 



O ' 1 INTRODUCTION 
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Narrow-line Seyfert 1 (NLS1) galaxies are a class of AGN 
characterized by a rather narrow width of the H/3 emission 
line (FWHM < 2000 km/s), and by the presence of a strong 
Fell bump, a soft X-ray excess and flux ratio [OIII]/H/3 <3 
jOsterbrock fc Poggdll985l ; |Poggell2000h . These sources are 
usually radio-quiet, although, in a few cases, they show 
a radio-loudness paramete r R (ratio between the 5 GHz 
and optical B flux densities, iKellermann et a l. 1989) greater 

20081 ). Such 
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than 1000 (|Komossa et al.l 
sources, dubbed radio-loud narrow-line Seyfert 1 galaxies 
(RL-NLS1), were thought to be inactive in 7-rays, although 
several a uthors speculated the occurrence of sim i larities with 
blaza r s (fzhou et al.ll20ol: iKomossa et al.l [200^ : IZhou et al ' 
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2JD10J). The important discovery of 7- ray emission from the 



RL-NLS1 source PMN J0948+0022 |Foschini et all l2010bl : 
Abd o et al . 2009a) confirmed these similarities, i.e. the pres- 
ence of a jet closely aligned to the line of sight as a source 
of Compton up-scattered 7-ray photons. PMN J0948+0022 
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(z=0.585, IWilliams et all I2002D is one of the strongest ra- 
dio source among the RL-NLS1; for its fast radio variabil- 
ity, source compactness, inverted radio spectrum and high 
brightness temperature, PMN J0948+0022 is one of the first 
RL-NL S1 for which the presence of a jet has been hypoth- 
esized (|Zhou et alJl2003l : lDoi et alj|2006h . The detection of 
7-ray emission definitively confirms the presence of a jet 
and allows to build a complete SED which closely resem- 
ble that of a typical blazar, with two broad peaks in the 
far IR and in the 7-ray range respectively. By modeling 
the SED with the synchrotro n and inverse-Compton model 
jGhisellini fc Tavecchic]|200g() it is possible to roughly esti- 
mate important parameters such as the black hole mass, the 
Eddington ratio, and the power carried by the jet. The re- 
sulting values, although model-dependent, are usually com- 
patible with estimates found independently in other studies. 
In the case of PMN J0948+0022 the black hole mass turns 
out to be ~1.5 x 10 s M^ (comp atible with the range of 
values found by IZhou et al.l 120031 using the empirical rela- 
tion between the FWHM of emission lines and the contin- 
uum luminosity), and the Eddington ratio which is ~40%. 
A 7-ray variability of a factor 2.2 rules out the possibil- 
ity t hat the 7-ray emiss ion is due to a starburst contribu- 
tion IjAbdo et al.ll2009bl ). Shortly after this discovery, three 
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more RL-NLS1 had been observed in 7-rays (|Abdo et al.l 
l2009d l : 1H 0323+342, PKS 1502+036 and PKS 2004-447, 
thus confirming that at least some RL-NLS1 are 7-ray emit- 
ting AGN. The SED of these sources has also been mod- 
eled with the synchrotron and inverse-Compton model, pro- 
viding black ho le mass estimates i n agreement with previ- 
ous studies (see lAbdo et al.ll2009c|, and ref erences therein) . 
1H 0323+342 (z=0.06. iMarcha et al.lll996l ) is considered to 
be a composite nucleus of a NLS1 (for its optical spec- 
trum properties typical of NLS1) and a blazar (for its flat 
radio spectrum, core compactness and X-rays variability, 
IZhou et al.l 120071 ). It is the only 7-ray detected RL-NLS1 
for which we have an _ffST/WFPC2 optical image of the 
host galaxy. The host shows a ring-like structure of 15.6 
kpc in diame ter and the entir e galaxy looks like a one- 
armed spiral (|Zhou et al.ll2007h . Another study, based on 
NOT data, suggests that 1H 0323+342 may be the remnant 
of a galaxy merger l|Anton et al.l I2OO8) . The most striking 
feature of this source is the high luminosity of the accre- 
tion disc which, according to S ED modeling, is ap proxi- 
mately 90% of Ed dington value (lAbdo et al.l 12009? ) . PKS 
1502+036 fz=0.41. ISnellen et al.ll2002l ) is the second most 
powerful source in terms of Eddington ratio (80 per cent). 
It carries a jet power, inferred from SED modelin g, sligthly 
lower than PMN J0948 +002 (lAbdo et al.l l2009d ). Finally, 
PKS 2004-447 (.2=0.24, iDrinkwater et all I199tF is an un- 
usual RL-NLS1. It is a strong radio emitter, with the highest 
radio-loudness parameter (R>6000) among the four 7-ray 
detected RL-NLS1. PKS 2004-447 has an unusually weak 
Fe II complex compared to other NLS1 and a rather steep 
(although varia ble) radio spec t ral in dex compared to the 
other RL-NLS1. lOshlack etail (|200ll ) suggests that it may 
be well classified as a low-luminosity compact steep spec- 
trum (CSS) radio quasar, a lso considering its low inferred 
black hole mass (~10 6 Mq) (|Abdo et al.ll2009cl ). Spectral fit 
of th e SED may requir e a thermal Comptonization compo- 
nent l|Gallo et al. 2006) or an external Compton component 
l|Abdo et al.ll2009d ). 

The detection of 7-rays from the four mentioned RL- 
NLS1 allows to identify a new class of 7-ray emitting AGN. 
The radio properties (namely temporal variability, flat spec- 
trum and high brightness temperature) together with the 
7-ray detection suggests the presence of a relativistic jet 
closely aligned to the line of sight. Low power, mildly rela- 
tivistic and poorly collimated radio jets have already been 
observed in a few spiral galaxies hosting Seyfert nuclei (e.g. 
iKeel et al]|2006l . and references therein), but in the case of 
PMN J0948+0022 and PKS 1502+036 the power carried by 
the jet, as inferred from SED modeling, is in the range of 
quasars, while in 1H 032 3+342 and PKS 20 04-447 is in the 
range of BL Lac objects ( Abdo et al. 2009c). Such powerful 
jets are observed on ly in blazars hosted in elliptical galax- 
ies l|Marscherl [2009] ) with black hole masses in the range 
10 8 - 10 9 Mq. By contrast black hole masses in 7-ray de- 
tected RL-NLS1 (10 6 - 10 8 Mq) are up to three orders of 
magnitude smaller than for blazars, thus suggesting that 
these sources emit at high Eddington ratios. Furthermore, 
it seems that the Mbh-o"* scaling relation does not apply to 
NLS1, mainly due t o their small FWHM of permitted lines. 
iDecarli et al.l f2008) showed that a reconciliation with the 
scaling relation is possible if the broad line region is assum ed 
to have a disc-like geometry (but see lMarconi et al"1l2008l ). 



While the 7-ray variability of PMN J 0948+0022 has al- 
ready been studied in lAbdo et all (|2009bl ) and, with a higher 
level of significance, during the outburst occurred in July 
2010 (|Foschini et alj l2010d ). it has never been studied for 
the remaining three sources. Aim of this paper is to ana- 
lyze the Fermi/LAT light curves of the afore-mentioned RL- 
NLS1 galaxies in order to put the 7-ray variability on a firm 
basis, and to find the minimum 7-ray variability timescale. 
Throughout the paper, we assume a ACDM cosmology with 
H = 71 km s" 1 Mpc" 1 , fi m = 0.27, Q A = 0.73. 



2 DATA ANALYSIS 

The data analysis has been performed following the stan- 
dard procedure described in the Fermi/LAT documentation. 
The Science Tools software version is 9.15.2 and the IRF is 
P6_V3_DIFFUSE. Data span the period from 2008 august 4 to 
2010 October 08 (~26 months). For each RL-NLS1 all class 

3 (diffuse) events coming from a zenith angle < 105° were 
extracted within an acceptance cone of radius 10° around 
the catalog source position. The unbinned likelihood analy- 
sis has been performed modeling the spectrum of each RL- 
NLS1 and of all nearby (< 10°) sources presen t in the LAT 
1-year point source catalog (|Abdo et al.ll2010bl ) with a power 
law in the range 0.1 - 100 GeV. We computed the integrated 
flux for each source analyzing all data over the entire period 
of 26 months. The extension of the region of interest to 15° 
around the catalog position results in non-significant varia- 
tions of the flux, n umber of counts and TS (Test Statistic, 
iMattox et ai1ll996l ), as expected since the PSF of LAT is 
at most 5° wide at 100 MeV. We also performed the same 
analysis using alternative spectral models. The use of a bro- 
ken power law instead of a simple power law model yields to 
a slight increase in the values of TS for PMN J0948+0022 
(from 993 to 1194) and 1H 0323+342 (from 86 to 115), and 
a decrease for PKS 1502+036 and PKS 2004-447. With a 
log-parabola, the values of TS remain essentially the same 
as for the case of a power law for PMN J0948+0022 (from 
993 to 999) and PKS 2004-447 (from 97.4 to 98.5), while 
they decrease for 1H 0323+342 and PKS 1502+036. In the 
following analysis we will therefore use the power law mod- 
els. Then we extracted light curves with time binning of 15 
days (Fig. Q] the time binning for PKS 2004-447 is 30 days) 
using a TS threshold o f 10 (TS>10, roughly equivalent to 
3cr, IMattox et ail [l996) and performed a chi-squared test 
against the null hypotesis of constant flux. If the detection 
was not significant (TS<10) we computed an upper limit to 
the flux by varying the source flux value (obtained through 
maximization of lik elihood) until TS reaches a value of 4 
l|Abdo et al.ll2010bl ). The resulting fluxes (denoted with ar- 
rows in the figure) corresponds to 2a upper limits. When 
TS<1 we didn't compute the upper limit since it would be 
overestimated. 

We further proceeded on the analysis of the light curves 
in order to compute the minimum e-folding timescale for 
each source. We extracted light curves with different time 
binnings starting from 30 days. When the detection is sig- 
nificant (TS>10) we re-run the analysis halving the time 
bin interval, down to a minimum of approximately 6 hours 
(roughly corresponding to four Fermi orbits, thus ensuring 
that each source is observed at least twice for each tem- 
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Figure 1. Upper panels: light curves of the four RL-NLS1 for detections with TS>10. Fluxes are given in units of 10 6 ph cm 2 s 1 
in the range 0.1 - 100 GeV. Vertical error bars correspond to lcr errors, while horizontal bars corresponds to the time binning (15 days 
for PMN J0948+0022, 1H 0323+342 and PKS 1502+036, 30 days for PKS 2004-447). Upper limits (TS<10) at 2a level are denoted 
by arrows. Middle panels: photon indices. Vertical error bars correspond to lcr errors. In both panels horizontal dashed lines are the 
integrated (over the entire period of 26 months) flux value and photon index respectively. Lower panels: TS values (plus symbols) and 
number of counts (triangle symbols, values on the right axis), the horizontal dotted line is the threshold (TS = 10). 



poral bin). Then, we considered all combinations of non- 
overlapping time bins with the following characteristics: (1) 
both time bins have a significant detection, TS>10; (2) the 
flux difference is greater than the greatest flux error involved 
at the 3<r level; (3) the count difference is significant at the 
3cr level, assuming a Poisson statistic; (4) the number of 
counts in each bin must be greater than 3. For such pairs of 
bins we computed the e-folding timescale as: 



Cj t j 



In Fi/Fj 



(1) 



where i and j are the indices of the involved time bins (i > 
j). The associated error (at 3a level) is computed through 
error propagation: 



An. 



{Atijf + T; 



AFj 

lJ F 



AFj 

' F, 



1/2 



(2) 

where Atij is half the width of the wider time bin involved, 
AFi and AFj are the 3cr error on fluxes. We compute the 
minimum e-folding variability timescale as r = min (fy), 
and associate the corresponding error At using Eq. [2] The 
reliability of Eq. [2] has been assessed by simulating 4 se- 
ries (respectively for ti, Fi, tj and Fj) of 10000 normally 
distributed values. The resulting values of r are normally 
distributed and the 3<r confidence intervals are correctly es- 
timated by At when the mean and standard deviation of 
the simulated data are used in Eq. [2] 



3 RESULTS 

Fig- Q] shows the light curves (upper panels) of the four RL- 
NLS1 over the entire period of 26 months. Vertical error bars 
are la errors on fluxes, while horizontal bars are equal to 
the time binning (15 days, except PKS 2004-447 for which is 
30 days) for all points with a significant detection (TS>10). 
Points without a significant detection are upper limits de- 
noted by arrows. Middle panels show the photon indices for 
points with significant detection, assuming a simple power 
law model (F(E) oc E F ) for each source in the range 0.1 - 
100 GeV. Horizontal dashed lines in both upper and middle 
panels are the integrated (over 26 months) flux and photon 
index respectively. Lower panels show the TS value (plus 
symbol) and number of counts (triangle symbol) for each 
bin with significant detection. The integrated luminosity 
and photon indices, as well as the \ 2 values and DOF in 
the hypotesis of constant flux, are reported in Table [T] The 
minimum e-folding variability timescale computed with Eq. 
[1] and [2] is reported in Table [1] The detailed view (zoom) 
on the time bins involved in the computation of the mini- 
mum e-folding variability timescale are shown in Fig. [2] To 
compute r, points with different time binnings are consid- 
ered and are denoted with different horizontal bar lengths 
in Fig. [2] (e.g. for PMN J0948+0022, one time binning is 
0.23 days and the other is 7.5 days. The bar corresponding 
to the shorter time bin is barely visible). Similarly, vertical 
bars refer to la errors in the flux. The square symbols iden- 
tify the points that fulfill the constraints described in Sect. 
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Figure 2. Detail view on the bins (denoted by squares) involved in the computation of the e- folding minimum variability timescale. 
Units and meaning of symbols are the same as in Fig. [l] (upper and lower panels) . Flux symbols have been changed to X for a clearer 
visibility. See also Table [2] 



Table 1. Data and results of the analysis on the four RL-NLS1 sources. Columns are: (1) name of the 
source; (2) redshift; (3) luminosity distance; (4) radio loudness; (5) integrated 7-ray luminosity (0.1 - 100 
GeV) over the entire period (26 months) with errors at lcr level; (6) photon index with errors at la level; 
(7) x 2 an d (8) DOF computed on the light curves of Fig. [T] in the null hypotesis of constant flux equal to 
the integrated flux; (9) minimum e-folding variability timescale with error at 3<r level. 



Source 


z 


D L 


R 




r 


x 2 


DOF 


T 






[Gpc] 




[10 45 erg s- 1 ] 








[days] 


PMN J0948+0022 


0.59 


3.40 


1000 


250.00 ± 13.04 


-2.851 ± 0.007 


528 


33 


3.3 ± 2.5 


1H 0323+342 


0.06 


0.27 


151 


0.16 ± 0.04 


-2.807 ± 0.010 


2575 


5 


17.7 ± 14.4 


my PKS 1502+036 


0.41 


2.20 


1549 


41.45 ± 4.10 


-2.708 ± 0.007 


385 


21 


12.0 ± 9.0 


PKS 2004-447 


0.24 


1.20 


6320 


3.85 ± 0.70 


-2.650 ± 0.006 


4032 


4 


28.4 ± 18.5 



Table 2. Quantities involved in the computation of the minimum e-folding variability timescale. 



Source 


C 






N / At 


TS 




[days] 


[days] 


[10~ 6 ph cm" 2 s" 1 ] 


[cts days -1 ] 




PMN J0948+0022 


703.20 


0.23 


2.578 ±0.584 


67.5 


38 




709.40 


7.50 


0.389 ±0.058 


7.0 


53 


1H 0323+342 


102.40 


0.94 


0.622 ±0.117 


21.5 


21 




135.70 


30.00 


0.095 ± 0.034 


2.6 


29 


PKS 1502+036 


198.00 


0.94 


0.674 ±0.058 


15.8 


16 




225.70 


30.00 


0.067 ± 0.025 


1.9 


11 


PKS 2004-447 


75.66 


30.00 


0.001 ±0.001 


0.1 


16 




214.40 


7.50 


0.197 ±0.006 


3.5 


11 



a MJD - 54682 
^ Time binning. 

c Flux in the range 0.1 - 100 GeV. 



© 2011 RAS, MNRAS 000, QHZ] 



Gamma-ray variability of radio-loud narrow-line Seyfert 1 galaxies 5 



[2] and provide the value of r. Notice, as an example, that 
the point at time ~706 (MJD - 54682, PMN J0948+0022) 
is excluded from the computation since its flux value is com- 
patible with nearby points at 3a level (rule (2) of Sect. [2]). 
The quantities involved in the calculation of r (time, time 
binning, photon flux, number of counts and TS values) are 
reported in Table [2] 



4 DISCUSSION 

A statistically significant variability over the entire period 
of 26 months is present for all sources. The significance of 
this variability on timescales of < 2 years is supported by 
the chi-squared test performed against the null hypotesis 
of constant flux (last two columns of Table [TJ. This rules 
out the possibility that the 7-ray emission is due to a star- 
burst activity. Thus, the data support the hypotesis that 
7-ray photons are associated to the presence of a jet. We 
cannot exclude the possibility of a starburst activity but its 
contribution would be negligible compared to the jet emis- 
sion, since the 7-ray luminosities (Table [T]) found in our RL- 
NLS1 are at least four order of magnitude greater than the 
arche typal starburst gala xy M82 hosting a quiescent black 
hole (|Gaffnev et al.ll 19931 ) and whose 7- ray luminosity in t he 
0.1 - 100 GeV range is ~10 40 erg s _1 l|Abdo et al.ll2010al) . 

The next step in our analysis has been to estimate the 
minimum timescale variability for each source. Among the 
many methods to measure a timescale variability we chose 
the e- folding timescale (Eq. [I}. The main advantage of this 
method is that it allows a computation of the timescale using 
just two flux measurements; it does not require any fitting or 
minimising procedure. Furthermore the resulting timescale 
is well defined and can be used to compare different sources. 
Finally, the error on the timescale is easily computed us- 
ing analytical error propagation (Eq. [2]). The underlying 
assumption is that the flux evolve according to an (either 
increasing or decreasing) exponential law. Although this as- 
sumption is not always justified (e.g. in the presence of flar- 
ing episodes) it is of common practice since we often do not 
know the actual law which drives the evolution of the flux. 
Furthermore, the requirement of just two flux measurements 
can be easily accomplished even for weak sources such as the 
ones analyzed in this work. The rules to select the flux mea- 
surements used in the computation of the e-folding timescale 
have been described in Sect. [5] rules (2) and (3) are needed 
to ensure that the flux in the two temporal bins are actually 
different, otherwise our assumption of exponentially varying 
flux cannot be justified. Rules (1) and (4) ensure that the 
flux measures and associated errors are reliable. The rela- 
tive error Ar/r of our measure of the minimum e-folding 
timescale depends on several factors. The main sources of 
uncertainties are the errors in the flux measurements AF 
and the width of the time bins At. The ratio AF/F can be 
made smaller by increasing the number of counts detected 
in a wider time bin (greater At). Viceversa, a narrower time 
bin results in a smaller number of counts and consequently 
a poor accuracy in the flux measurement. The best achiev- 
able accuracy on r (that is the smaller value of Ar/r) is 
thus determined by the trade-off between the narrowness of 
the time bin and the accuracy of flux measurements. Notice 
however that the narrowness of the time bin is limited by 



the requirement that the source is detected with high sig- 
nificance, TS> 10 (rule (1) in Sect. [2]). Since the shortest 
width of the time bins scales linearly with the inverse of 
the photon flux, we may use narrower time bins (and de- 
tect shorter timescales) on brighter sources. The accuracy 
of r is improved if we compute the e-folding timescale us- 
ing non-contiguous time bins, since the exponential law is 
better constrained by distant points rather than closer ones. 
On the other hand, if the bins are contiguous, the error 
Ar may be significantly greater than r itself, i.e. the mea- 
sure would be useless. The issue related to the computation 
of r with non-contiguous time bins is that we are deliber- 
ately ignoring what lies between the two time bins, even 
if we have a significant measure. An example is given in 
Fig. [2] for PMN J0948+0022 (upper-left panel) where the 
minimum e-folding timescale is computed ignoring the flux 
measurement at time ~ 706 (MJD - 54682). The reason 
is that the flux measure has a rather big error bar, which 
does not allow us to consider it statistically incompatible 
at the 3<r level with neighbouring flux measures (rule (2) 
in Sect. [2]). This example suggests that the intrinsic vari- 
ability timescale, i.e. a measure of "how fast" the source is 
able to change its flux significantly, may be shorter than 
our measure of r, although we cannot assess it on a firm 
statistical basis. This is a consequence of the fact that we 
are assuming a particular function (the exponential) and 
that we are using only two flux measurements instead of 
all available data. If our detector were able to monitor the 
flux evolution continuously, we could have employed more 
sophisticated methods, such as the Fourier analysis, in or- 
der to estimate the intrinsic minimum timescale variability. 
Actually, our light curves do not allow a continuous mon- 
itoring of the sources with time binning of 30 days, and 
we are thus forced to rely on the e-folding timescale. As a 
consequence, our measure of r should be considered as an 
upper limit for the intrinsic minimum timescale variability. 
The sources display an e-folding minimum timescale vari- 
ability in the range of 3 - 30 days. In particular, variability 
of the order of -days has been found for PMN J0948+0022, 
and of the order of -tens of days for 1H 0323+342, PKS 
1502+036 and PKS 2004-447. The minimum measured vari- 
ability timescales in blazars can be as low as 200 seconds 
at very high energies (E>200 G eV) and 800 seconds at X- 
rays (see I Aharonian et al.ll2007l . and references therein). At 
Fermi/LAT energies variability on sc ales of few hours has 
been detected in several blazars (e.g. Foschini et alj 120081 ; 
iTavecchio et all 120101 : iFoschini et alJl2010al ).' In particular, 



the bright blazars 3C 454.3 and PKS 1510-089 showed vari- 
ability on scales of 3-6 hours with flux — 10~ 5 ph cm -2 s _1 
(|Tavecchio et all 120101 ). that is 1 - 2 orders of magnitude 
greater than the flux of the sources analyzed here. As dis- 
cussed above, we expect to measure longer timescales on 
weaker sources, as a consequence of the longer time inte- 
gration required to significantly detect the source. Our min- 
imum timescale estimates are indeed approximately 1-2 
orders of magnitude longer compared to those of the above- 
mentioned blazars. Therefore, we cannot exclude the possi- 
bility that variability in RL-NLS1 is as fast as the one ob- 
served in some of the most luminous blazars. The detection 
of shorter timescales on our sources is challenging due to 
their weakness. Ground-based Cherenkov telescopes may be 
more suited to search for shorter timescale variabilities, if the 
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spectrum of RL-NLS1 extend to very high energies. In prin- 
ciple variability can also be measured using the flux upper 
limits and the resulting minimum timescales would be be- 
low the results quoted in Table [Q (i.e. for PMN J0948+0022 
we would obtain r ~ 1 day). Upper limit of fluxes are how- 
ever less reliable than direct flux measurement, since they 
depends on the contribution of nearby sources and of diffuse 
background. 

The minimum e-folding timescale variability allow to 
estimate an upper limit for the size of the emitting region: 
flwob < 5cr/(l + z) -0.5 - 6 <5 x x I0 17 cm, where Si = 5/10 
is the relativistic Doppler factor and z is the redshift. With 
a jet aperture 8-i = 9/0.1 rad the distance at which most 
of the kinetic energy of the jet is dissipated and the 7-rays 
observed are produced (the so-called dissipation region), is 
^?diss < 0.2 - 2 Si/0-i pc. The photon indices are rather 
steep (r < —2, with the only exception of the point at time 
~75 of PKS 2004-447) , thus the inverse Compton peak lies 
below 100 MeV. In the framework of the blazar sequence 
l|Fossati et alj|l998l ) the 7-ray emitting RL-NLS1 are there- 
fore located in the region relevant to quasars. The photon 
index shows also some significant variations, e.g. at time 
~750 (MJD - 54682) the photon index of 1H 0323+342 un- 
derwent a 6<r variation from V = - (2.46 ±0.06) to F = - 
(3.2 ± 0.1) in 15 days (TS = 15 and 24, counts = 43 and 74 
respectively). The corresponding photon flux is larger when 
the spectrum is steeper, so that the overall 7-ray luminosity 
does not change significantly. Another striking case is given 
at time ~550 (MJD - 54682) of PKS 1502+036 for which the 
photon index shows a 16cr variation from V — -(1.97 ± 0.06) 
to T = -(4.7 ± 0.2) in 15 days, although the TS values are 
much lower (10 and 11 respectively). Also in this case the 
overall 7-ray luminosity does not change significantly. This 
behaviour is in contradiction with the harder-when-brigther 
phenomenon observed in the X-ray spectra o f several HBL 
i|Brinkmann et al.l [20051 ; ISembav et alll2002ft and at 7-ray 
energies of PMN J0948+0022 (|Foschini et aJj|2010d ). This 
phenomenon is usually ascribed to the upshift of the inverse- 
Compton peak as the source brightens. The weakness of the 
sources prevents us from building a time-resolved detailed 
spectra and from drawing any conclusion about the eventual 
shift towards lower energies of the inverse-Compton peak in 
1H 0323+342 and PKS 1502+036. 



5 CONCLUSIONS 

In this work, we report the discovery of a statistically sig- 
nificant variability in the 7-ray light curves of the four RL- 
NLS1 detected with Fermi/LAT , with minimum variability 
timescales in the range 3-30 days. This excludes a potential 
starburst origin of the 7-ray emission, and supports the hy- 
pothesis of the presence of a jet closely aligned to the line of 
sight. A hint for photon index variations on timescales ~tens 
of days is also found in the data. Variability appears to be 
a feature common to the four first 7-ray detected RL-NLS1. 
The minimum timescales found in the Fermi/LAT energy 
range, appropriately scaled with the flux, are comparable to 
those found in the most luminous blazars. Thus, it is not 
possible to exclude variability as fast as that observed in 
bl azars. This study goes in the same direction of the finding 
by iFoschini et alj 2010d who reported compelling evidence 



of similarities in the SED shape of PMN J0948+0022 (also 
analyzed here) and the archetypal blazar 3C 273. We are 
confident that in the near future more RL-NLS1 will be 
identified as 7-ray emitters among the many unidentified 
7-ray sources observed by Fermi/LAT . 



REFERENCES 

Abdo A. A., et al, 2009a, ApJ, 699, 976 
Abdo A. A., et al., 2009b, ApJ, 707, 727 
Abdo A. A., et al., 2009c, ApJ, 707, L142 
Abdo A. A., et al., 2010a, ApJ, 709, L152 
Abdo A. A., et al., 2010b, ApJS, 188, 405 
Aharonian F., et al., 2007, ApJ, 664, L71 
Anton S., Browne I. W. A., Marcha M. J., 2008, A&A, 490, 
583 

Brinkmann W., Papadakis I. E., Raeth C, Mimica P., 

Haberl F., 2005, A&A, 443, 397 
Decarli R., et al., 2008, MNRAS, 386, L15 
Doi A., et al., 2006, PASJ, 58, 829 

Drinkwater M. J., Webster R. L., Francis P. J., Condon 
J. J., Ellison S. L., Jauncey D. L., Lovell J., Peterson 
B. A., Savage A., 1997, MNRAS, 284, 85 
Foschini L., Longo F., Iafrate G. et al., 2008, ATel 1784 
Foschini L., et al., 2009, Advances in Space Research, 43, 
889 

Foschini L., et al., 2010a, MNRAS, 408, 448 

Foschini L. for the Fermi/LAT Collaboration, Ghisellini G., 
Maraschi L., Tavecchio F., Angelakis E., 2010b, in "Accre- 
tion and ejection in AGN: a global view" , eds L. Maraschi, 
G. Ghisellini, R. Delia Ceca & F. Tavecchio, Como (Italy), 
22-26 June 2009, ASP Conference Series 427, p. 243 

Foschini L., et al., 2010c, MNRAS submitted, 
arXiv: 1010.4434 

Fossati G., et al., 1998, MNRAS, 299, 433 

Gaffney N. I., Lester D. F., Telesco C. M., 1993, ApJ, 407, 
L57 

Gallo L. C, et al., 2006, MNRAS, 370, 245 
Ghisellini G., Tavecchio F., 2009, MNRAS, 397, 985 
Gu M., Chen Y., 2010, AJ, 139, 2612 
Keel W. C, et al., 2006, AJ, 132, 2233 
Kellermann K., et al., 1989, AJ, 98, 1195 
Komossa S., et al., 2006, AJ, 132, 531 

Marcha M. J. M., Browne I. W. A., Impey C. D., Smith 

P. S., 1996, MNRAS, 281, 425 
Marconi A., et al., 2008, ApJ, 678, 693 
Marscher A., 2009, In: "The Jet Paradigm - From Micro- 
quasars to Quasars", ed T. Belloni, Lect. Notes Phys. 794, 

in press, (arXiv: 0909. 2576) 
Mattox J. R., et al., 1996, ApJ, 461, 396 
Oshlack A. Y. K. N., Webster R. L., Whiting M. T., 2001, 

ApJ, 558, 578 
Osterbrock D. E., Pogge R. W., 1985, ApJ, 297, 166 
Pogge R. W., 2000, New Astronomy Review, 44, 381 
Sembay S., Edelson R., Markowitz A., Griffiths R. G., 

Turner M. J. L., 2002, ApJ, 574, 634 
Snellen I. A. G., McMahon R. G., Hook I. M., Browne 

I. W. A., 2002, MNRAS, 329, 700 
Tavecchio F., et al., 2010, MNRAS, 405, L94 
Williams R. J., Pogge R. W., Mathur S., 2002, AJ, 124, 

3042 



© 2011 RAS, MNRAS 000, [THT] 



Gamma-ray variability of radio-loud narrow-line Seyfert 1 galaxies 



Yuan W., et al., 2008, ApJ, 685, 801 
Zhou H., et al., 2003, ApJ, 584, 147 
Zhou H., et al., 2007, ApJ, 658, L13 

This paper has been typeset from a TffjX/ WF^X. file prepared 
by the author. 



© 2011 RAS, MNRAS 000, [THT] 



